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We study the low energy states of Kondo alloys as function of the magnetic impurity concentration per site, 
x, and the conduction electron average site occupation, n c . Using two complementary approaches, the mean- 
field coherent potential approximation and the strong coupling limit, we identify and characterize two different 
Fermi liquid regimes. We propose that both regimes are separated by a Lifshitz transition at x = n c . Indeed, 
we predict a discontinuity of the number of quasiparticles which are enclosed in the Fermi surface. This feature 
could provide a scenario for the non-Fermi liquid properties that were recently observed in Kondo alloy systems 
around x = n c . 
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Kondo alloys are realized in broad families of strongly cor- 
related materials where magnetic quantum impurities are em- 
bedded randomly in a crystalline system with conduction elec- 
trons [1]. The Kondo problem has been intensively studied 
for almost half a century [1], and the single impurity model 
was solved exactly by various methods [2-6]. In the 70's, 
Nozieres has adapted Landau's approach [7, 8], showing that 
a single Kondo impurity can be described universally at low 
energy as a local Fermi liquid tJ. Multiple-impurity models 
have later been introduced in order to describe systems with 
a high concentration of Kondo ions. It appeared that, unlike 
the single impurity models, the multiple-impurity models are 
not belonging to a single universality class, but can rather give 
rise to a large variety of ground states depending on the lattice 
structure, the electronic occupation, and on the strength of the 
Kondo coupling which can be tuned with pressure 
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Putting aside the relevant issues of ordering, we concen- 
trate here on the Heavy Fermi liquid ground state that is ob- 
served at low temperature in lots of dense Kondo systems Jl2ll . 
The electronic exhaustion problem suggested bv 

Nozieres remained a long standing issue for years ll 15142311 
and was finally understood |2ll [24]: the Kondo tempera- 
ture Tk that characterizes the crossover to the low temper- 
ature Kondo screening regime, and the Fermi liquid energy 
scale T co h that characterizes the Fermi liquid ground state 
are indeed the same energy scale. It was shown that the ra- 
tio Tk I T co h still depends on the electronic filling and lattice 
structure |2ll [25l |26|]. But this ratio was proven not to de- 
pend on the Kondo coupling as it was wrongly thought ini- 
tially El EH. 

A simple picture of the periodic dense Kondo lattice 
model had been provided within the strong Kondo coupling 
limit fl^[l^l27l - [29ll . Even if this limit does not correspond 
to the experimental reality, its qualitative validity is supported 
by poor man's scaling analysis til QOfl, which shows that 
the low temperature Kondo physics of the single impurity 



model is renormalized to a strong coupling fixed point. In 
the infinite coupling limit, the ground state of a Kondo lat- 
tice made of N sites with Ns = N spins 1/2 is charac- 
terized by the formation of N c local Kondo singlets, where 
N c denotes the number of conduction electrons. The coher- 
ent Fermi liquid is recovered from a perturbation expansion 
at lowest order in hopping /coupling. The corresponding 
fermionic quasiparticles correspond to the bachelor Kondo 
spins, whose double occupancy is forbidden, giving rise to 
strong correlation effects. This strong coupling description 
also provides the correct number of quasiparticles: N — N c 
free spins, which, according to a particle-hole transformation, 
give 2N - (N — N c ) = N c + N quasiparticles (the factor 
2 comes from the spin-degeneracy). This result is in perfect 
agreement with Luttinger theorem Bill , which predicts an en- 
largement of the Fermi surface due to the contribution of the 
Kondo impurities 11141 12111 . 

Our analysis starts here: on the one side the universal single 
impurity Kondo model can be described as a local Fermi liq- 
uid. This picture describes well Kondo systems with very di- 
lute impurities. On the other side, many realizations of dense 
Kondo systems with a periodic lattice of magnetic ions have 
a Fermi liquid ground state. Can these two Fermi liquids be 
connected continuously to each other? 

We consider a Kondo alloy model, defined by the following 
Hamiltonian 
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where c^J denotes annihilation (creation) operator of a con- 
duction electron with spin a =t, I on site i of a periodic lat- 
tice that contains TV sites. Uj denotes intersite hopping en- 
ergy, and Jk is a local Kondo antiferromagnetic interaction 
between local quantum spin 1/2 denoted Si and the local den- 
sity of spin of conduction electrons, <Ji. These Ns = xN 
Kondo spins, with concentration x < 1, are located on a sub- 
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part K of the periodic lattice. We choose a frozen configu- 
ration for K— sites, whose positions are distributed randomly 
without spatial correlation. The number of conduction elec- 
trons is fixed to N c = n c N, and we assume for simplicity that 
the electronic filling satisfies n c < 1. 

First, the Kondo alloy model is studied using the mean- 
field approximation for the Kondo interaction llil [32[ [33|] 
and considering the random position of K— sites within a 
CPA-DMFT method l25ll . i.e., a matrix version of Dynami- 
cal Mean-Field Theory (DMFT) (33j] which is equivalent 
to the matrix version of the Coherent Potential Approxima- 
tion (CPA) [37|] . Kondo alloys models had been studied 
earlier using CPA methods, but focusing mostly on the high 
concentration of impurities Jlsl-lii. The matrix CPA-DMFT 
formalism which is presented in details in Ref [25] treats sep- 
arately Kondo and non-Kondo sites. For the sake of simplicity 
the numerical calculations were performed assuming the con- 
duction electrons move on a Bethe lattice, characterized by an 
elliptic non-interacting local density of states (d.o.s.), as de- 
picted on figure[T] The local d.o.s. of the conduction electrons 
on a K— site is analyzed as functions of x, n c , and the temper- 
ature, T. At high T (left column on figure [B it corresponds 
to the non-interacting conduction band, and for T & Tk (sec- 
ond column) a Kondo resonance forms around the Fermi level 
uo = 0. The position of the resonance depends on n c , but 
its shape at T « Tk does not depend on x. This is fully 
consistent with the interpretation that Tk characterizes the in- 
coherent local Kondo singlet formation. At lower temperature 
(right column), when the Fermi liquid regime sets in, the local 
d.o.s. can exhibit two qualitatively different features, depend- 
ing on the impurity concentration: for x < n c (figures Eh 
andHt), it remains roughly similar to the non-interacting one, 
with just a small deformation around the Fermi level, due to 
the Kondo resonance. But for x > n c (figures and Hi) 
the Kondo resonance is split by an hybridization gap. Note 
that a reduction of the density of state in a dense Kondo sys- 
tem had already been suggested in 13 8. kill. This splitting can 
be explained by the standard two-band picture which is com- 
monly used for Kondo lattices. A systematic analysis of the 
local d.o.s. at low T for different impurity doping and elec- 
tronic filling clearly shows that the crossover between dense 
(with a gap) and dilute (without gap) regimes occurs when 
x = n c . The dense system remains Fermi liquid since the 
gap occurs at uj > 0, i.e., above the Fermi level. However, 
this mean-field result suggests that the Kondo-interacting and 
non-interacting Fermi liquids can be smoothly connected to 
each other only for x < n c . At larger impurity concentrations 
(or, equivalently, at lower electronic filling), a singularity may 
disconnect the two Fermi liquid regimes, related to this gap 
opening. 

Studying the strong coupling limit, Jk ^> \Uj\, provides 
a clear evidence of the analytical disconnection between the 
coherent and the local Fermi liquids which are expected in the 
dense and the dilute regimes. Hereafter we generalize the ex- 
act strong coupling treatment of Ref. 11511 to the Kondo alloy 
model. 
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FIG. 1. Local d.o.s. on a Kondo site for electronic filling n c = 0.20 
(a and b) and n c = 0.60 (c and d), and impurity concentrations a) 
x = 0.10, b)x = 0.30, c)x = 0.50, and d) x = 0.70. Left column: 
elliptic d.o.s. obtained at T ^> Tk. Second column: formation of 
the Kondo resonance, at T « Tk. Right column: comparison be- 
tween the dilute local and dense coherent Fermi liquid regimes, with 
respectively x < n c (a and c), and x > n c (b and d). Insets: hy- 
bridization gaps above the Fermi level, which characterize the coher- 
ent dense regime. Parameters: Jk = 1.5 and non interacting elec- 
tronic bandwidth W = 4. d.o.s. for T <C Tk have been computed 
with T ^T K /1000, with T K = 0.099 and 0.18 for n c = 0.20 and 
0.60 respectively. 
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FIG. 2. Schematic picture of the ground state of the Kondo alloy 
model in the large Jk limit. Here, we chose N = 12 and iV c = 6 
(green simple arrows). Sites with Kondo spins are depicted by red 
double arrows. Black points refer to non-Kondo sites. Kondo singlets 
are marked by a solid circle, a) Dilute case with Ns = 2 < iV c : 
quasiparticles are identified to the extra iV c — N s electrons moving 
freely on the non-Kondo sites, b) Dense case with Ns = 8 > iV c : 
quasiparticles are identified to the Ns — N c unscreened Kondo spins 
(doted circles) moving on the Kondo sites. 



We start with the extreme limit of tij = 0. Depending 
on the ratio N c /N$ = n c /x, we can distinguish two differ- 
ent situations, as illustrated on figure [2 In the dilute regime 
Ns < N c the ground state is characterized by Ns local 
Kondo singlets on each Kondo- site. Ns electrons are thus 
localized on the K— sites which effectively become "forbid- 
den sites" for the remaining N c — Ns electrons. The latter 
can freely occupy the left N — Ns non-Kondo sites. This de- 
generate ground state has a finite entropy per site S dilute = 
(1 — x) m(l — x) — (n c — x) ln(n c — x) — (l — n c ) m(l — n c ). 
In the dense regime Ns > N c the ground state is charac- 
terized by the formation of 7V C Kondo singlets. The non- 
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Kondo sites are empty, but a degeneracy results from the var- 
ious energetically equivalent possible positions of the Kondo 
screened impurities on K— sites. The resulting entropy per 
site is Sdense = xln(x) - n c m(n c ) - (x - n c ) \n(x - n c ). 
Both S dilute and Sdense are analytical functions of x and n c , 
but they cannot be analytically connected to each other. In- 
deed, introducing the relative density S = x - n c , we find 
that both vanish like S « -\5\ In \S\ in the limit 5^0. The 
presence of an absolute value in this expression is clearly iden- 
tified as a singularity separating the dense and dilute regimes. 
Furthermore, the fact that the infinite coupling entropy van- 
ishes at Ns = N c suggests that this point will not correspond 
to a Fermi liquid even when considering perturbations in elec- 
tronic hopping t. 

Then, we consider the lowest order in perturbation expan- 
sion, which is linear in t. Fermi liquids characterize both di- 
lute and dense regimes, but a singular point separates the local 
Fermi Liquid from the coherent Fermi liquid. 

Dilute regime Ns < N c : the first perturbative correction 
generates an effective model where the fermionic quasiparti- 
cles are N^ lute = N c — Ns free electrons moving with hop- 
ping energy Uj on the lattice made of the N — Ns non-Kondo 
sites (see figure |2] a). Of course this is valid only if 1 — x is 
bigger than the percolation threshold x p . The only correlation 
effect in this effective model results from the frozen random 
depletion of the Kondo- sites. This picture recovers the early 
approaches of Nozieres to the single impurity model, which 
was mapped onto a problem of light conduction electrons scat- 
tering on the Kondo singlet. Here, light conduction electrons 
scatter on the K— sites which are effectively excluded. From 
CPA approach, we find that the site depletion leads to an ef- 
fective conduction band for the N^ lute quasiparticles, with a 
renormalized bandwidth corresponding to a rescaled hopping 
term y/l — x. The local d.o.s is mostly renormalized by the 
site depletion, and its total spectral weight is reduced by a fac- 
tor 1 — x. The missing weight, x, is transfered to a higher en- 
ergy peak, which corresponds to excitations where one of the 
K- sites has or 2 conduction electrons. Note that the occupa- 
tion Nqp lute = N c — Ns corresponds to a filling of (n c — x)/2 
quasiparticles per site and per spin component; this quantity 
is smaller than 1 — x since we assumed an electronic filling 
n c < 1 (a particle-hole transformation would generalize this 
to n c > 1). 

Dense regime Ns > N c : performing perturbation ex- 
pansion at lowest order we obtain an effective model where 
fermionic quasiparticles correspond to the unscreened Kondo 
impurities. The mapping is very similar to the one described 
for a periodic Kondo lattice Ill5ll . Here, we assume that the 
K— sites concentration exceeds the percolation threshold, i.e., 
that x > x p . The unscreened impurities move on a depleted 
lattice made of the K— sites. There are Ns — N c unscreened 
spins moving on Ns sites, and since there motion is hole-like, 
the number of quasiparticles is N^ nse = 2N S - (N s -N c ) = 
Ns + N c . This suggests that the Luttinger theorem applies in 
the dense regime, both electrons and Kondo spins contribut- 
ing to the Fermi surface. Here, the correlations have two ori- 
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TABLE I. Characteristics of the dense and dilute regimes in the 
strong coupling limit. 



gins: the random depletion, which also characterizes the di- 
lute regime, and a supplementary infinite repulsion which pre- 
vents from unphysical double impurity occupation on a same 
site. Since the moving holes are Kondo singlets, the lattice 
hopping is also renormalized by a factor 1 / 2 as explained in 
Ref. H13TI . An extra renormalization of the hopping results 
from the lattice site depletion. Invoking a CPA approach, the 
intersite hopping for the correlated quasiparticles is Ujy/x/2. 
The local d.o.s. is renormalized in two steps: first, the site de- 
pletion reduces the bandwidth by a factor y/x, and the spec- 
tral weight by a factor x. The missing spectral weight here 
corresponds to excitations moving an electron from a Kondo 
singlet to a non Kondo site (with or 2 conduction electrons). 
Then, the peak at lowest energy is split again due to the ef- 
fective Hubbard repulsion. We may interpret the upper for- 
bidden Hubbard band as corresponding to excitations of local 
Kondo singlets into local triplets. This singlet- triplet excita- 
tion mode remains gaped (or pseudo-gaped) in the opposite 
limit, Jk <C t, as indicated by mean-field calculations (see 
figure[T]). Note that the occupation N^ nse = Ns + N c corre- 
sponds to a filling of (x + n c )/2 quasiparticles per site and per 
spin component; this quantity is smaller than x/2 in the dense 
regime, and it is associated to a chemical potential inside the 
lowest Hubbard band. This mapping to an effective Hubbard 
model is exact in the strong Kondo coupling limit, but it is in- 
teresting to observe that the numerical results obtained from 
the mean-field approximation with weak coupling parameters 
also provide a splitting of the d.o.s. with a gap opening in the 
dense regime (see figure [T]). 

The characteristics of the strong coupling effective models 
obtained in the dense and dilute regimes are compared in ta- 
ble U The impossibility of connecting analytically these two 
families of models thus survives to finite t corrections. We 
identify N c = Ns to a critical point with the following singu- 
larities: the number of quasiparticles is discontinuous, as well 
as the number of effective sites, the effective hopping and the 
Hubbard repulsion, which changes from to oo. Also the 
Fermi surface collapses in a discontinuous way. For N c = N, 
this singular point corresponds to a Kondo insulator, and is 
physically accessible only from the dilute side. For N c < N, 
the singularity of this point remains, and it is accessible from 
both dilute and dense sides. The critical point is not associated 
to a symmetry breaking, but to a singular change of the Fermi 
surface known as a Lifshitz transition (LT) 114211 . 
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Number of Kondo spins (Ns) 



FIG. 3. Evolution of the number of quasiparticles as a function of 
the number of Kondo spins Ns for a fixed occupation of conduction 
electrons N c . A transition is expected at Ns — N c , when x p < 1/2 
(see table HU). 
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TABLE II. Situations where Kondo disorder and Lifshitz transitions 
can be observed. 



The evolution of the number of quasiparticles is shown in 
figure [3] as a function of the number of Kondo spins in the 
strong coupling limit. The depicted situation corresponds 
to a Kondo alloy which is realized by substitution of non- 
Kondo (example La) with Kondo (example Ce) atoms, with 
fixed light electron structure and filling. In the dense impurity 
regime, the Kondo spins contribute to the volume of the Fermi 
surface, which is consistent with the general idea of Luttinger 
theorem. The formation of this large Fermi surface requires a 
sufficiently high amount of Kondo impurities, in order to form 
a coherent heavy Fermi liquid. In the dilute regime, Luttinger 
theorem is expected to break down for Ns < N c . There, 
Kondo impurities are too diluted for contributing to the Fermi 
surface. Instead, their main effect is rather to capture con- 
duction electrons for the formation of localized Kondo sin- 
glets. Only the resulting N c — Ns electrons contribute to the 
Fermi surface, and scatter on the localized Kondo singlets. 
The discontinuity of N qp is related to a Lifshitz transition: the 
Fermi surface is not defined at the critical point N c = Ns. 
This is true for N c = Ns = 1 because the system is in this 
case a Kondo insulator, but we predict this situation to occur 
more generally for N c = Ns not necessarily equal to 1. Ap- 
proached from the dilute side, this critical point corresponds to 
a continuous collapse of the Fermi surface. Approached from 
the dense side, the collapse of the Fermi surface is discontinu- 
ous. Of course, this analysis does not consider possible issues 
of cluster formations and we also assumed percolation of the 
depleted effective lattices. Various possible cases are summa- 
rized in table III An intermediate Kondo disorder (KD) regime 
with cluster formation is expected for lattices with big perco- 
lation threshold, x p > 1/2 (which is usual in two-dimensional 
systems with only nearest neighbor hopping). The most real- 
istic cases are three-dimensional systems with low percolation 



threshold x p < 1/2 and intermediate electronic filling. Two 
situations can thus occur: the Lifshitz transition should be ob- 
served for x p < n c < 1 — x p , whilst a Kondo disorder phase 
is expected for other values of n c . 

Here, we did not analyze the Kondo disorder regime, where 
inhomogeneities may control the physical properties as dis- 
cussed in Ref. L43]. Even without considering these possi- 
ble other effects, the predicted collapse of the Fermi surface 
suggests that the physical properties of a Kondo alloy in the 
vicinity of the critical point Ns = N c should be non Fermi 
liquid (NFL). This mechanism may provide a possible sce- 
nario for the origin of the NFL properties observed recently 
in the heavy fermion compounds Ce ;E Lai_ ;E Ni2Ge2 114411 and 
Ce^Lai-^Ptln [45]. A difference (transition or crossover) be- 
tween dilute Kondo system and a dense heavy fermion had 
also been suggested by other theoretical methods, including 
finite size calculations [29], quantum Monte Carlo simula- 
tions [46, 47] and Gutzwiller approximation [48ll. We pro- 
pose to analyze more systematically Kondo alloys in order 
to find other examples of NFL behavior of this type. These 
studies should be completed by measurements which ana- 
lyze the Fermi surface (Angleresolved Photoemission Spec- 
troscopy, quantum oscillations), in order to test the predicted 
validity and violation of Luttinger theorem, respectively in the 
dense and dilute regimes. One difficulty that may be faced 
in some compounds would come from a complex electronic 
band structure. Usually, several conduction bands are in- 
volved, and the criteria N c = Ns might be not always reach- 
able. This critical point could also have signatures in other 
experimental realizations of Kondo alloys, in bulk materials 
or in artificial atom systems like quantum dot arrays or optical 
lattices. 
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